Intrinsic rotation of toroidally confined magnetohydrodynamics 
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The spatiotemporal self-organization of viscoresistive magnetohydrodynamics (MHD) in a toroi- 
dal geometry is studied. Curl- free toroidal magnetic and electric fields are imposed. It is observed 
in our simulations that a flow is generated, which evolves from dominant ly poloidal to toroidal 
when the Lundquist numbers are increased. It is shown that this toroidal organization of the flow 
is consistent with the tendency of the velocity field to align with the magnetic field. Up-down 
asymmetry of the geometry causes the generation of a non-zero toroidal angular momentum. 
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a. Introduction. The magnetic confinement of fu- 
sion plasmas is strongly influenced by turbulent fluctu- 
ations. These fluctuations degrade the quality of the 
confinement and thereby reduce the performance of the 
fusion reactor. It was discovered three decades ago [1] 
that, under certain circumstances, the turbulent activity 
is reduced, leading to a better confinement. Still today 
the understanding of this low-to-high-confinement (LH) 
transition is far from complete. There is however strong 
evidence that large toroidal velocities of the plasma are 
a feature that is either at the origin, or a consequence of, 
the mechanism that is responsible for the LH-transition 
0, E|. Large toroidal velocities, of the order of several 
kilometers per second, are observed even in the absence of 
external momentum input. Several mechanisms are put 
forward to explain the toroidal rotation, mostly based 
on the turbulent transport of toroidal momentum gen- 
erated at the tokamak edge (e.g. in references [4-6]). In 
this letter we present a mechanism which seems to be 
generic, since it is observed even in one of the coarsest 
descriptions of a fusion plasma: visco-resistive magneto- 
hydrodynamics (MHD). 



b. A MHD description of fusion plasmas. In the 
MHD description, the plasma is described as a charge- 
neutral conducting fluid. MHD, despite its low level of 
complexity compared to kinetic descriptions or two-fluid 
descriptions, gives already rise to a wealth of intricate 
phenomena and its analytical treatment is only possible 
in some simplified cases, either in the absence of velocity 
fields [3, [8j or in the absence of nonlinear interactions [9| . 
We will come back to these analytical approaches, but 
before that, we present the equations that we consider. 
These are the dimensionless incompressible viscoresistive 
MHD equations for the velocity field u and for the mag- 
netic field £?, in Alfvenic units [10j . 



— - M _1 V 2 n 
at 



dB 
~dt 
E 
V • u 



- V ( P + -u 2 ) + u x :uj + j x B, 



-V x E, 

— [uxB] , 



= 0, 



V-B = 0, 



(1) 

(2) 

(3) 
(4) 




Ro I R 




Asymmetric 'D' 



Circular 



FIG. 1: Cross-sections of the toroidal geometries considered 
in the present work. The toroidal angle is labelled and the 
poloidal one ip. 



with the current density j = V x £?, the vorticity 
cj = V x n, the pressure P and the electric field 
E. These equations are non-dimensionalized using the 
toroidal Alfven speed Ca — Bo / y/p/no as typical veloc- 
ity, with Bq the reference toroidal magnetic field at the 
center of the torus (R = Rq), p the density and //q the 
magnetic constant. The reference length L (see Fig. [T]) 
is the diameter of the cross section for the circular case 
and is the minor diameter for the asymmetric 'D' shape 
(L = 1.88 for both geometries). The dynamics are then 
governed by the initial and boundary conditions of the 
problem, and two dimensionless quantities: the viscous 
Lundquist number (M) and the Lundquist number (S) 
defined as 
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with A the magnetic diffusivity and v the kinematic vis- 
cosity. The ratio of these two quantities is the magnetic 
Prandtl number P r = v/\, which we have chosen unity 
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in the present study, thereby reducing the number of free 
parameters, which characterize the magnetofluid, to one, 
the viscous Lundquist number. Previous investigations 
indicate that it is the geometric mean of the viscosity 
and the magnetic diffusivity which determines the dy- 
namics [Tl],u3- I n setting the Prandtl number to one, a 
change in the Lundquist numbers, M or 5, is equivalent 
to a change in the Hartmann number. 

Let us now go back to the analytical description of 
visco-resistive MHD. In the static case in which u = 0, 
equation (pQ) reduces to an equilibrium 

VP = j x B. (6) 

In a cylindrical geometry this equilibrium can be achieved 
by various magnetic configurations such as the z-pinch or 
the #-pinch [13[ . In toroidal geometry it is problematic to 
obtain such an equilibrium, as we will now explain. We 
consider the case in which the driving toroidal electric 
field is curl- free within the plasma, over times of interest, 
such that Eq ~ 1/R. Further do we assume the toroidal 
magnetic field to obey the same scaling, which follows 
from the integration of Ampere's law on a toroidal loop. 
In the simplest case, we choose a space- uniform electrical 
conductivity such that the toroidal current induced by 
the electric field is also given by the same dependence, so 
that the externally imposed magnetic field and toroidal, 
laminar, voltage-driven current density are given by, 

B (R) = Bo^-e 9 , MR) = Jo^e 6 . (7) 

Computing the Lorentz force resulting from these 
toroidal fields, taking into account the poloidal mag- 
netic field induced by Jo, results in a force-field which 
is not curl- free [l4|. Since the curl of the pressure gra- 
dient is necessarily zero, the equilibrium described by 
(j6j) becomes impossible and additional terms of equation 
(P) need to be taken into account to balance the equa- 
tion. Since all other terms in (pp) are proportional to (or 
quadratic in) the velocity, the resulting state must be 
dynamic. That is, a toroidal plasma, described by visco- 
resistive MHD, confined by curl- free toroidal electric and 
magnetic fields, necessarily moves. 

It is true that the rationale described above depends 
on the choice of the electric conductivity, which was as- 
sumed to be uniform. It was however shown [l5|, [l6| 
that to satisfy © in a torus, very unusual profiles of 
the electrical conductivity must be assumed. We omit 
these rather unphysical cases and focus on the dynamical 
plasma behaviour which results for the simplest, uniform, 
conductivity profile. 

It follows from the foregoing that it is necessary to 
take into account all other terms in the MHD equa- 
tions, and analytical treatment becomes impossible un- 
less symmetries are assumed. To study the full dynam- 
ics we are obliged to solve numerically the system and 
this is what is done in the present investigation. Such 
fully three-dimensional non-stationary simulations, tak- 
ing into account all relevant time and space-scales, are 




FIG. 2: Streamlines colored by the value of the toroidal veloc- 
ity, uq for M = 15 (top) and M = 150 (bottom) in the geom- 
etry with circular cross-section. Only a part of the toroidal 
domain is shown. 



computationally demanding and only quite recently have 
the necessary resources and numerical methods become 
available to do such simulations. Equations (pE]) are 
discretized with a Fourier pseudo-spectral method on a 
Cartesian grid. To impose the boundary conditions we 
use the volume-penalization technique, a method of the 
immersed boundary type, which we consider a good com- 
promise between the ease of implementation, flexibility in 
geometry, and the numerical cost of the simulation. Re- 
sults for two dimensional MHD can be found in reference 
[I?}- We recently extended this method to study the 
three-dimensional visco-resistive MHD equations [l8j], 
and in the present communication we present the results 
of three-dimensional simulations in two toroidal geome- 
tries. 

c. Results of numerical simulations. Details of the 
numerical method are given in reference [18| . Simulations 
are carried out on a cubic domain of size 2tt consisting 
of 256 3 grid-points for the highest values of M. The 
initial condition for the simulations is zero velocity, and 
no-slip velocity boundary conditions are imposed. We 
consider the boundaries of the fluid domain as perfectly 
conducting and coated with an infinitely thin layer of in- 
sulator. Thereby the normal component at the wall of the 
magnetic and current density fields vanishes. We impose 
toroidal magnetic and current density fields given by eq. 
0. the Biot-Savart law is used to determine the poloidal 
magnetic field induced by the toroidal current Jq(R). All 
the simulations presented in this communication are per- 
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FIG. 3: The ratio of the mean-square toroidal velocity to the 
total mean-square velocity (uq) / (|^ 2 |) as a function of M. In 
the inset we show the average over the domain of the absolute 
value of the cosine of the angle between the velocity field and 
magnetic field. 



formed with B = 0.8 and Jo = 0.3. This corresponds, 
for both geometries, to a pinch ratio ~ 0.16, defined 
as the ratio between the wall-averaged poloidal and the 
volume- averaged toroidal magnetic field (Q = B^/ (Be)). 
The only parameter that we vary is the Lundquist num- 
ber M . The simulations are time-dependent and they are 
stopped when a dynamical steady state is reached. 

The results in Fig. [2] show the presence of a poloidal 
flow, a pair of counter-rotating vortices in the r — cp plane. 
For small M the dynamics are dominantly poloidal, as is 
expected. Indeed, in the limit of vanishing nonlinear- 
ity, Bates et al. [9| showed analytically that the steady 
state solution is a pair of poloidally rotating vortices, 
aligned with the toroidal direction. For non-zero nonlin- 
earity, i.e. , by increasing M, the vortices start moving in 
the toroidal direction, both in opposite direction. Their 
toroidal velocity increases with the Lundquist number M 
in the two considered geometries. The three-dimensional 
velocity streamlines show a substantial change of topol- 
ogy from dominantly poloidal to dominantly toroidal flow 
(see Fig. [2j bottom). This is quantified in Fig. [3j where 
we observe that the principal direction of the flow mo- 
tion is toroidal if M is raised beyond ~ 40. The square 
toroidal velocity saturates for increasing M at a value 
of ~ 86% of the total square velocity. This toroidal or- 
ganization of the flow is consistent with the tendency of 
the velocity-field to align with the magnetic field, as is 
illustrated in the inset of Figure [3j where we compute the 
average (over the toroidal domain) of the absolute value 
of the cosine of the angle between the velocity and mag- 
netic field. This quantity is equal to one if the velocity 
and magnetic field are perfectly aligned or anti-aligned. 
It is shown that the trend towards a toroidal velocity fol- 
lows exactly the same M dependence as the alignment. 

This tendency towards a dominant toroidal flow is sim- 



ilar for the torus with the asymmetric cross-section as is 
shown in Fig. [3j while the alignment is even more pro- 
nounced. In both geometries, the generated velocity field 
contains non- negligible fluctuations. The quantities u' 
and B' denote the fluctuations around the azimuthally 
averaged instantaneous velocity and magnetic field, re- 
spectively. At M = 3008 in the D-shaped geometry, 
Kms/^rms = 6.7 • HT 2 , B' r7ns /B rrns = 1.5 • 10" 3 and 
u rms /B rms = 3.2 • 10 -3 . The rms values correspond 
here to volume averages over the toroidal domain. A 
detailed investigation of the spatial distribution of these 
fluctuations and its dependence on M will be presented 
elsewhere. 

A fundamental difference is observed between the flows 
that are generated in the two geometries. The volume- 
averaged toroidal angular momentum is defined by 

(L e ) = ^ J^RuedV. (8) 

For the torus with circular cross-section, the time-average 
of this quantity is zero to a good computational approxi- 
mation, due to the up-down symmetry of the observed 
flow. However, for the torus with asymmetric cross- 
section this is not the case. For low M, a poloidal pair 
of counter-rotating vortices appears (see Figure 2] for az- 
imuthally averaged flow visualizations, and toroidal ve- 
locity profiles along a vertical cut) as for the circular 
cross-section. Similarly, if the viscous Lundquist number 
is increased, an important toroidal flow develops. Unlike 
the symmetric case, there is a breaking of the symme- 
try in the flow and the part of the flow moving in the 
negative direction (the red zone), becomes larger at the 
expense of the part of the flow which moves in the positive 
toroidal direction (blue zone). This symmetry breaking, 
illustrated in Fig. 2J leads to the development of a net 
toroidal flow. The toroidal angular momentum becomes 
hereby non-zero (see Fig. [5j). Its normalized value in- 
creases significantly with the viscous Lundquist number. 
The observed influence of up-down symmetry is consis- 
tent with axi-symmetric time-independent computations 
[lOj and is also observed in gyrokinetic simulations and 
experiments [19|, [20|. 

It is presently not clear if the velocity profile observed 
in our simulations will change qualitatively when M is in- 
creased further and a transition to another flow-topology 
cannot be excluded. Also have we not yet investigated 
the influence of the magnetic Prandtl number. It is at 
this point perhaps important to say that we do not know 
what the viscosity should be to approximate the dynam- 
ics of experiments. However, the fact that this feature 
is observed in fully resolved simulations of the visco- 
resistive MHD equations is a result of major importance, 
since it shows how intrinsic toroidal rotation is present in 
one of the coarsest global descriptions of a fusion plasma, 
without invoking arguments on charge non-neutrality or 
kinetic theory. 

d. Conclusion. We want to summarize the results 
that we obtained: considering curl-free toroidal electric 
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FIG. 5: Normalized toroidal angular momentum 
\{Lq)\ / Le rrns as a function of M observed in the tori 
with asymmetric and symmetric cross-section, respectively. 

and magnetic fields and constant transport coefficients, 
visco-resistive magnetofluids spontaneously generate ve- 



locity fields. This velocity field aligns (or anti-aligns) 
with the magnetic field, thereby generating a toroidal 
component. This is a nonlinear effect which becomes 
negligible in the limit of small Lundquist number. Fur- 
thermore, toroidal angular momentum is created, if the 
up-down symmetry of the torus is broken. 

By its simplification, both in terms of the used model- 
equations, as well as in terms of the parameter range 
chosen for the properties of the conducting fluid, our in- 
vestigation should be considered academic rather than 
directly applicable to the detailed description of exist- 
ing machines. At the same time, since MHD does give 
a rough description of laboratory plasmas, the mecha- 
nism that we have described should be present, at least 
qualitatively, in existing devices. The observed magne- 
tohydrodynamic self-organization thereby seems to be of 
major importance for the magnetically confined fusion 
community. 

This work was supported by the contract SiCoMHD 
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IDRIS, project 22206. 
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